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1. Introduction 



The muon g — 2 has served over many decades as an important test of the Standard Model 
(SM). It is also sensitive to contributions from New Physics slightly above the electroweak scale. 
In fact, for several years now a discrepancy of more than three standard deviations has existed 
between the SM prediction and the experimental value, see the recent reviews Refs. [@, |, [| on the 
muon g — 2. The main error in the theoretical SM prediction comes from hadronic contributions, 
i.e. hadronic vacuum polarization and hadronic light-by-light (had. LbyL) scattering. Whereas the 
hadronic vacuum polarization contribution can be related to the cross section e + e~ — > hadrons, 
no direct experimental information is available for had. LbyL scattering. One therefore has to rely 
on hadronic models to describe the strongly interacting, nonperturbative dynamics at the relevant 
scales from the muon mass up to about 2 GeV. This leads to large uncertainties, see Refs. [[|, [|] 
for recent reviews on had. LbyL scattering, largely based on the original works [||, |7|, |8j, |9p. 

Essentially, these models describe the interactions of hadrons with photons, usually with the 
help of some form factors. One can reduce this model dependence and the corresponding un- 
certainties by relating the form factors at low energies to results from chiral perturbation theory 
(ChPT) @ and at high energies (short distances) to the operator product expansion (OPE) [jnj]. 
In this way, one connects the form factors to the underlying theory of QCD. In particular, this has 
been done in Refs. [H [7, 12, || ||] for the numerically dominant contribution from the exchange of 
light pseudoscalars % , r\ , r\' . 

The pseudoscalar-exchange contributions to had. LbyL scattering are given by the diagrams 
shown in Fig. [l[ 



% 7i°, r|,T|' 





Figure 1: The pseudoscalar-exchange contributions to had. LbyL scattering. The shaded blobs represent the 



form factor ^ P s*Y"f wnere PS = 7T ,tj,tj', 7Z° 



It was pointed out recently in Ref. [g], that one should use fully off-shell form factors for the 
evaluation of the LbyL scattering contribution. This seems to have been overlooked in the recent 
literature, in particular, in Refs. [12, ||, |9[ The on-shell form factors as used in Refs. []8| [l2|] 



actually violate four-momentum conservation at the external vertex, as observed already in Ref. [g]. 
The exchange of the lightest state 71° yields the largest contribution and therefore warrants 



special attention. In these proceedings, based on the results obtained in Ref. [|13|], we present a 
new QCD short-distance constraint on the off-shell pion-photon-photon form factor J^o*^ at 
the external vertex by relating it to the quark condensate magnetic susceptibility of QCD. We then 



evaluate this contribution in the framework of large-Afc QCD [14], using a form factor which fulfills 
this new and other relevant short-distance constraints. 
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2. On-shell versus off-shell form factors 



For the pion, the key object which enters the diagrams in Fig. [I] is the off-shell form factor 
^^^.^.{{qi + <?2 ) 2 , ^ i , ^1 ) which can be defined via the QCD Green's function (VVP) [|6|, fj[ |l3|] 



rf 4 xrf 4 ye ! ' (91 " +?2 - y) (0|r{j M (x) iv (y)P 3 (0)}|0) 
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up to small mixing effects with the states f] and r\' and neglecting exchanges of heavier states like 
7iP , n° , — Here (x) is the light quark part of the electromagnetic current and P 3 = T//7/5 4- 
The corresponding contribution to the muon g — 2 may be worked out with the result [||] 
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where the external photon has now zero four-momentum. See Ref. [Jsj] for the expressions for 7}. 

Instead of the representation in Eq. (2.2), Refs. [JTJ, |8|] considered on-shell form factors which 
would yield the so called pion-pole contribution, e.g. for the term involving T 2 , one would write [Q] 

x ^ Jl o rr (ml,(q l +q 2 ) 2 ,0). (2.3) 



^n a 7* 7* ( m fl; ) 01 1 02 ) 



Although pole dominance might be expected to give a reasonable approximation, it is not correct as 
it was used in those references, as stressed in Refs. ^ g]. The point is that the form factor sitting at 
the external photon vertex in the pole approximation J^w^m 2 ., (q\ +q 2 ) 2 ,0) for (171 +<?2) 2 7^ w| 
violates four-momentum conservation, since the momentum of the external (soft) photon vanishes. 
The latter requires ^ n Q*y*y{{q\ +^2) 2 , (01 +02) 2 ,O). Ref. [^] then proposed to use instead 

.^■fy-{m 2 n ,q\,q\) x J^o^ra 2 ,m 2 ,0) . (2.4) 

Note that putting the pion on-shell at the external vertex automatically leads to a constant form 
factor, given by the Wess-Zumino-Witten (WZW) term jjTJ]. However, this prescription does not 



yield the pion-exchange contribution with off-shell form factors, which we calculate with Eq. (2.2). 



Strictly speaking, the identification of the pion-exchange contribution is only possible, if the 
pion is on-shell. If one is off the mass shell of the exchanged particle, it is not possible to separate 
different contributions to the g — 2, unless one uses some particular model where elementary pions 
can propagate. In this sense, only the pion-pole contribution with on-shell form factors can be 
defined, at least in principle, in a model-independent way. On the other hand, the pion-pole contri- 
bution is only a part of the full result, since in general, e.g. using some resonance Lagrangian, the 
form factors will enter the calculation with off-shell momenta. In this respect, we view our eval- 
uation as being a part of a full calculation of had. LbyL scattering using a resonance Lagrangian 
whose coefficients are tuned in such a way as to systematically reproduce the relevant QCD short- 



distance constraints, along the lines of the resonance chiral theory developed in Ref. [16] 
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3. A new short-distance constraint on the off-shell pion-photon-photon form factor 



The form factor J^o-y^^gi + <72) 2 ><7i><72) defined in Eq. (2.1) is determined by nonpertur- 



bative physics of QCD and cannot (yet) be calculated from first principles. Therefore, various 
hadronic models have been used in the literature. At low energies, the form factor is normalized 
by the decay amplitude, j# (71° — > 77) = e 2 ^ r o rr (m|,0,0). To a good approximation, all hadronic 
models thus have to satisfy the constraint j^o^fm 2 ,0,0) = —Nc/(12n F n ). 

For an on-shell pion, there is also experimental data available for one on-shell and one off-shell 
photon, from the process e + e~ — > e + e~K° '. Several experiments Jl9| ] thereby fairly well confirm 
the Brodsky-Lepage j j20| ] behavior for large Euclidean momentum lim^p^^ ^ r ^o r * r (m|, — Q 2 ,0) ~ 
-2F n /Q 2 and any model should reproduce this behavior, maybe with a different prefactor. 

Apart from these experimental constraints, any consistent hadronic model for the off-shell 
form factor ^^^((qx +^) 2 , ^1,^2) should match at large momentum with short-distance con- 
straints from QCD that can be calculated using the OPE. In Ref. [|^] the short-distance properties 



for the three-point function (VVP) in Eq. ( |2.l| ) in the chiral limit and assuming octet symmetry 
have been worked out in detail. Two limits are of interest. In the first case, the two momenta 
become simultaneously large, which describes the situation where the space-time arguments of all 
three operators tend towards the same point at the same rate. The second situation corresponds 
to the case where the relative distance between only two of the three operators in (VVP) becomes 
small. When the space-time arguments of the two vector currents in (VVP) approach each other, 
the leading term in the OPE leads to the Green's function (AP). The explicit results for both these 



cases can be found in Refs. [22, 13] 



The new short-distance constraint on the off-shell form factor at the external vertex in had. 
LbyL scattering arises when the space-time argument of one of the vector currents in (VVP) ap- 
proaches the argument of the pseudoscalar density. This leads to the two-point function (VT) of 
the vector current and the antisymmetric tensor density 

^(nvrWOO = / d 4 xe^(0\T{V^x)(W<y p aYV)(0)}\0), ° P o = ±[Y P ,Yo]- (3-1) 

Conservation of the vector current and invariance under parity then give (TL\j)^ pa (p) = (p p r\^ a — 
/?c7 T J,up)nvT(/? 2 ). In this way one obtains the relation (up to corrections of order a s ) p2| , [13| ] 



hm M ^o Vr 4(A^+^ 2 ) 2 ,(Ag 1 ) 2 ,^) = -^^^n V T(^) + ^({j • (3.2) 

In particular, at the external vertex in LbyL scattering in Eq. (^2|), the limit ^2 — ► is relevant. 
As pointed out in Ref. [23], the value of Yhji(p 2 ) at zero momentum is related to the quark 



condensate magnetic susceptibility % in QCD in the presence of a constant external electromagnetic 



field, introduced in Ref. [£4|]: (0\qa^ v q\0) F = ee q x(Y¥)oFfiv, with e u = 2/3 and e rf = -1/3. 



1 We note that in our work Jl3| and in Refs. [Q, ^, |§| ^ simply F K = 92.4 MeV is used, without any error attached. 
Maybe this could be an additional source of uncertainty in a^ yh,7t , in particular in view of the new value r(7T° — » 77) = 
(7.82 ± 0.23) eV presented in Ref. @; see also the discussion in Ref. ^ and references therein. 

2 Note, however, that a recent measurement of the form factor by the BABAR collaboration [ pT| ] at momentum 
transfers Q 2 between 4 GeV 2 and 40 GeV 2 does not show such a falloff . We will come back to this issue in Section 0. 
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With our definition of ITvt in Eq. (3.1) one obtains the relation nyr(O) = —((WY)o/2)X ( see a l so 



Ref. [25]) and the new short-distance constraint at the external vertex can be written as [|13l 



lim J^o 



7i" y 7 



((^,) 2 ,(A„ 



(3.3) 



Note that there is no falloff in this limit, unless % vanishes. 

Unfortunately there is no agreement in the literature what the actual value of % should be. In 
comparing different results one has to keep in mind that % actually depends on the renormalization 
scale /J. In Ref. [Q] the estimate = 0.5 GeV) = -(8.16l 2 ^) GeV 2 was given in a QCD sum 
rule evaluation of nucleon magnetic moments. A similar value x = — Nc/(4n 2 F 2 ) = —8.9 GeV 2 
], probably again for a low scale ji. ~ 0.5 GeV as argued in Ref. [26]. 



was obtained in Ref. 

On the other hand, saturating the leading short-distance behavior of the two-point function 
with one multiplet of lowest-lying vector mesons (LMD) [28, 23, 22] leads to the estimate 
^lmd _ —2/ My = —3.3 GeV -2 @]. Again, it is not obvious at which scale this relation holds, 
it might be at /i = My. This LMD estimate was soon afterwards improved by taking into account 
higher resonance states (p', p") in the framework of QCD sum rules, with the results % (0.5 GeV) = 
-(5.7 ±0.6) GeV" 2 [El and #(1 GeV) = -(4.4 ±0.4) GeV 



A more recent analysis H3 



yields, however, a smaller absolute value #(1 GeV) = —(3.15 ±0.30) GeV -2 , close to the original 
LMD estimate. 3 For a quantitative comparison of all these estimates for % we would have to run 
them to a common scale, for instance, 1 GeV or 2 GeV, which can obviously not be done within 
perturbation theory starting from such low scales as /I = 0.5 GeV. 



4. New evaluation of the pseudoscalar-exchange contribution in large-Afc QCD 

In the spirit of the minimal hadronic Ansatz j j33| ] for Green's functions in large-Afc QCD, an 
off-shell form factor JP n a* r -, r »((q\ + qij 1 ,q\,q\) has been constructed in Ref. [p2|]. It contains the 
two lightest multiplets of vector resonances, the p and the p' (LMD+V), and fulfills all the OPE 
constraints discussed earlier: 

fir q\ql{q\ + ql + ql)+PH{q 2 yql,ql) f4n 

3 (q 2 -Ml) (q 2 -M 2 2 ) (q 2 -M 2 ) {q\ -M 2 2 ) ' 

h {q\ + q2) 2 + h 2 q\ q\ + % {q\ + ql)ql + h * qt 

+h 5 (q 2 l +q 2 2 )+h 6 q 2 3 +h 1 , q\ = {q x + q 2 ) 2 . (4.2) 

Below we reevaluate the pion-exchange contribution using off-shell LMD+V form factors at 
both vertices. The constants hi in the Ansatz for J^ ¥^ v in Eq. ( |Q| ) are determined as follows. 
The normalization with the pion decay amplitude 71° — > 77 yields h~i = —NcM v M Vi J (4n 2 F 2 ) — 
h 6 ml - h 4 ml = - 14.83 GeV 6 - /j 6 m 2 - /j 4 «4 . where we used M Vl =M p = 775 .49 MeV and M Vl = 
M p i = 1.465 GeV. The Brodsky-Lepage behavior can be reproduced by choosing h\ = GeV 2 . 

3 After the publication of our paper Ref. Jl3|], two new estimates for % appeared, both based on the analysis of the 
zero-modes of the Dirac operator. Ref. | |3l| ] presents an analytical approach which yields ^(1 GeV) = —3.52 GeV~ 2 
with an estimated error of 30 — 50%. A quenched lattice calculation [^] for Nc = 2 gives a very small absolute value 
X = —1.547(6) GeV -2 . No scale dependence is given, the lattice spacing corresponds to 2 GeV. 



^^y^^ \q\q\^q 1 i) — 
P v H {qUlql) = 
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In Ref. [Q a fit to the CLEO data [0] for the on-shell form factor ^y?, +V (m| ; -Q 2 ,0) was 
performed, with the result h$ = (6.93 ±0.26) GeV 4 — h^m\. The constant I12 can be obtained from 
higher-twist corrections in the OPE with the result I12 = —10.63 GeV 2 

Within the LMD+V framework, the vector-tensor two-point function reads [ ^2[ [l3| ] 



n 



LMD+V 
VT 



(P 2 



(W>o 



(p2 



•M, 



CVT 



V 2 J 



As shown in Ref. [E3], the OPE constraint from Eq. fO) for ^ L ¥°+ v leads to the relation 



h[ +^3 +^4 = 2cvT- 



(4.3) 



(4.4) 



The LMD estimate £ LMD = -2/Af£ = -3.3 GeV" 2 is close to x(H = 1 GeV) = -(3.15 ± 
0.30) GeV" 2 obtained in Ref. [|(]] using QCD sum rules with several vector resonances p,p', 
and p" . Assuming that the LMD/LMD+V framework is self-consistent, we will therefore take 
X = ( — 3.3 ±1.1) GeV" 2 in our numerical evaluation, with a typical large-A^ uncertainty of about 
30%. We will vary hj, in the range ±10 GeV 2 and determine J14 from Eq. ( |4.4| ) and vice versa. 

The coefficient he in the LMD+V Ansatz is undetermined as well. It enters at order p 6 in 
the low-energy expansion of (VVP) in one combination of low-energy constants from the chiral 
Lagrangian of odd intrinsic parity, A^md+v = -F^h 6 /(SM^M^ 2 ) @]. The LMD ansatz with only 



-0.26 (10" 4 /F 2 ) ph. If the 



one multiplet of vector resonances yields A^^ +q \2 = —F%/(8My) 
LMD/LMD+V framework is self-consistent, the change in these estimates, while going from LMD 
to LMD+V, should not be too big. Since the size of this low-energy constant seems to be small 
compared to another combination of low-energy constants which enters at order p 6 , we allow for a 
100% uncertainty of A^P + )2 and get the range h 6 = (5 ± 5) GeV 4 , see Ref. [||] for details. 

The results for a^ >ylj ' 7t for some selected values of , h% and h^, varied in the ranges discussed 
above, for % = -3.3 GeV" 2 , hi = GeV 2 , h 2 = -10.63 GeV 2 and h 5 = 6.93 GeV 4 - h 3 ml are 
collected in Table [T], see Refs. [|T3|,[3|] for details on the numerics. 





h 6 = GeV 4 


he = 5 GeV 4 


h e = 10 GeV 4 


h 3 = -10 GeV 2 


68.4 


74.1 


80.2 


h 3 = GeV 2 


66.4 


71.9 


77.8 


h 3 = 10 GeV 2 


64.4 


69.7 


75.4 


h 4 = -10 GeV 2 


65.3 


70.7 


76.4 


h 4 = GeV 2 


67.3 


72.8 


78.8 


h 4 = 10 GeV 2 


69.2 


75.0 


81.2 



Table 1: Results for a^ hyUK ° x 10 11 obtained with the off-shell LMD+V form factor for % = -3.3 GeV" 2 
and the given values for hj,,h 4 and he- When varying hi (upper half of the table), the parameter h 4 is fixed 



by the constraint in Eq. (4.4). In the lower half the procedure is reversed. 



/y Lb vL " 7T^ 1 1 

Varying % by ±1.1 GeV" changes the result for by ±2.1 x 10" at most. The 

uncertainty in he affects the result by up to ±6.4 x 10" 11 . The variation of aj^ yL ' K with h 3 [with 
h 4 determined from the constraint inEq. ( |4.4| ) or vice versa] is much smaller, at most ±2.5 x 10" 11 . 
In the absence of more information on the values of the constants hi,,h 4 and he, we take the average 
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of the results obtained with h(, = 5 GeV 4 for h?, = GeV 2 and for h\ = GeV 2 as our central value: 
^LbyL.jr _ 3 x irj- 11 . Adding all uncertainties from the variations of %, hi, (or h^), h$ and 
linearly to cover the full range of values obtained with our scan of parameters, we get []l3|, ||] 

4 byL;7t ° = (72 ± 12) x HT 11 . (4.5) 

This value replaces the result obtained in Ref. [j|] with on-shell LMD+V form factors at both 
vertices. We think the 16% error should fairly well describe the inherent model uncertainty using 
the off-shell LMD+V form factor. In order to facilitate updates of our result in case some of the 



parameters hj in the LMD+V Ansatz in Eq. ( ]4. 1[ ) will be known more precisely, we have given in 



the Appendix of Ref. [13] a parametrization of a^ yL ' K for arbitrary coefficients hi. 4 

As far as the contribution to from the exchanges of the other light pseudoscalars v\ and r]' is 
concerned, it is not so straightforward to apply the above analysis within the LMD+V framework 



to these resonances. In particular, the short-distance analysis in Ref. [ |22| ] was performed in the 
chiral limit and assumed octet symmetry. We therefore resort to a simplified approach which was 
also adopted in other works !§, [7], ||, and take a simple VMD form factor normalized to the 
experimental decay width T(PS — > 77). In this way we obtain the results aj? yUn = 14.5 x 10" 11 
and a^ byL ' ?7 = 12.5 x 10~ n , which update the values given in Ref. [§]. Adding up the contributions 
from all the light pseudoscalar exchanges, we obtain the estimate [|l3|, ^ 

a LbyL;PS = (99±16)>< ^-11^ (4 6) 

where we have assumed a 16% error, as inferred above for the pion-exchange contribution. 5 



5. Discussion and conclusions 

We would like to stress that although our result for the pion-exchange contribution is not 
too far from the value a^ ylj,7t pole = (76.5 ±6.7) x 10 -11 given in Ref. [0], this is pure co- 
incidence. We have used off-shell LMD+V form factors at both vertices, whereas the au- 
thors of Ref. ^ evaluated the pion-pole contribution using the on-shell LMD+V form factor 
^■ n Qy>y*{q\,q\) = J^iyy^m^g 2 ,^) at the internal vertex and a constant WZW form factor at 
the external vertex, see for instance Eq. (18) in Ref. Since only the pion-pole contribution 
is considered in Ref. [Qj, their short-distance constraint cannot be applied to our approach either. 
However, our ansatz for the pion-exchange contribution agrees qualitatively with the short-distance 
behavior of the quark-loop derived in Ref. [Q], see the discussion in Refs. Jl3| , g]. 

4 A fit of the on-shell LMD+V form factor to the recent BABAR data Q yields h\ = (-0. 17 ± 0.02) GeV 2 and 
/J5 = (6.51 ±0.20) GeV 4 — h-jm^ with ^ 2 /dof = 15.0/15 = 1.0. In this way we would get the new average value 
^LbyL.jr _ g x 10~ n , i.e. the result is essentially unchanged fromEq. 

Applying the same procedure to the electron, we get a^ byUK = (2.98 ±0.34) x 10~ 14 |^. This number super- 
sedes the value given in Ref. Note that the naive rescaling ag byL ' Jt (rescaled) = (m e /m il ) 1 a^ yh ' n = 1.7 x 10~ 14 
yields a value which is almost a factor of 2 too small. Our estimates for the other pseudoscalars contributions us- 
ing VMD form factors at both vertices are ag byL ' T ' = 0.49 x 10~ 14 and a^ hyh '^ = 0.39 x 10~ 14 . Therefore we get 
^LbyL.PS = (3.9±0.5) x 10 — 14 , where the relative error of about 12% is again taken over from the pion-exchange contri- 
bution. Assuming that the pseudoscalar contribution yields the bulk of the result of the total had. LbyL scattering correc- 
tion, we obtain ag byL ' had = (3.9± 1.3) x 10~ 14 , with a conservative error of about 30%, see Ref. [g]. This value was later 
confirmed in the published version of Ref. |B| where a leading logs estimate yielded ag byL ' had = (3.5 ± 1.0) x 10~ 14 . 
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Our results for the pion and the sum of all pseudoscalars are about 20% larger than the values 
in Refs. [^J 0] which used other hadronic models. An evaluation of the pion-exchange contribution 
using an off-shell form factor based on a nonlocal chiral quark model yielded a^ yL ' K = (65 ±2) x 
10~ n [34]. In that model, off-shell effects of the pion always lead to a strong damping in the form 
factor and the result is therefore smaller than the pion-pole contribution obtained in Ref. [gp. In 
our model, there are some corners of the parameter space where the result is larger than the pion- 
pole contribution, for instance, we get a maximal value of a^ byL jr = 83.3 x 10~ n in the scanned 
region. Very recently, a value of a^ yL ' PS = 107 x 10~ n with an estimated error of at most 30% 



was obtained in Ref. [35] within an AdS/QCD approach. 

Combining our result for the pseudoscalars with the evaluation of the axial- vector contribution 
in Ref. ^ and the results from Ref. [||] for the other contributions, we obtain the estimate JT3| , ||] 

fl LbyL;had = (116±4())xl(r ll (51) 

for the total had. LbyL scattering contribution to the anomalous magnetic moment of the muon. To 
be conservative, we have added all the errors linearly, as has become customary in recent years. 
In the very recent review ^ the central values of some of the individual contributions to had. 
LbyL scattering are adjusted and some errors are enlarged to cover the results obtained by various 
groups which used different models. The errors are finally added in quadrature to yield the estimate 
^LbyL.had _ ±26) x 10~ n . Note that the dressed light quark loops are not included as a 
separate contribution in Ref. |@]. They are assumed to be already covered by using the short- 
distance constraint from Ref. [Q] on the pseudoscalar-pole contribution. Certainly, more work on 
the had. LbyL scattering contribution is needed to fully control all the uncertainties. 
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